Electron-positron pair creation by Coulomb and laser fields in the tunneling regime 
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Electron-positron pair creation due to combined nuclear Coulomb and strong laser fields is inves- 
tigated for the tunneling regime. The energy spectra and angular distributions of the pair are found 
analytically for the first time. The energy spectrum for each lepton exhibits a sharp maximum 
located well above the threshold for any polarization of the laser field. The angular distributions 
of leptons depend on the polarization: for the linear polarization both leptons move predominantly 
along the laser beam direction; for the circular polarization leptons are emitted in a thin-walled 
cone centered on the laser beam. The spectral and angular distributions found are governed by the 
intensity and frequency of the field, and the frequency independent total pair creation rates comply 
with the previously known results. A new method of calculation - the vtcinal approximation - which 
uses the fact that the pair production takes place in the close vicinity of the nucleus, is suggested. 

PACS numbers: 12.20.Ds, 32.80.Wr, 42.50.Hz 
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I. INTRODUCTION 

We consider the electron-positron pair creation due to 
a strong laser field and the Coulomb field of a nucleus. 
This combination of fields may create electron-positron 
pairs via a multiphoton process, sometimes called the 
nonlinear Bethe-Heitler process [ij. Such a scheme of 
electron-positron pair creation has received attention 
over the last decade due to advances in laser technol- 
ogy, which may permit it to be tested in an x-ray free 
electron laser (XFEL) experiment • 

The multiphoton pair creation by combined strong 
laser and Coulomb fields was first treated by Yakovlev Q . 
He considered total and partial cross sections of the pair 
creation process in the case of a circularly polarized laser 
field. Mittleman [ij later calculated the cross-section and 
pair creation rate for the low intensity linear polariza- 
tion case. More recently, Refs. d, B 0| have considered 
the closely related problem of multiphoton pair creation 
due to a ultra-relativistic nucleus colliding with a strong, 
circularly polarized laser field. The energy spectra and 
angular distributions of the produced leptons were calcu- 
lated numerically over a range of Lorentz factors of the 
nucleus. Ref. [7] also treated the linear polarization case 
and the scenario in which the pair is created in a bound 
state of the nucleus, and Ref. [SJ considered pair creation 
in the vicinity of heavy ions. Ref. @ has derived the to- 
tal creation rate due to a circularly polarized laser field 
incident on a nucleus and also calculated numerically the 
lepton energy spectrum for strong fields. 

One approach widely used in the literature treats the 
interaction of leptons with the laser field exactly by use of 
the Volkov states 0, [ll| . It is analogous to the Keldysh 
approximation used for the problem of multiphoton ion- 



isation of atoms or ions [13, M^M, [IE [H, [13, E [H 



Alternatively, Milstein et al. [20| have recently derived 
pair creation rates and cross-sections for different cases 
of polarizations and different parameters of the laser field 
by application of the optical theorem and use of the po- 
larization operator of a photon in a laser field. 

Overall, the total rates of the pair creation in the prob- 
lem considered have been studied in detail by different 
methods, which all result in clear, simple analytical ex- 
pressions. However, the more detailed characteristics of 
the problem, namely the spectral and angular distribu- 
tions of the leptons, have not been examined with the 
same level of detail. Meanwhile, considering future ex- 
perimental studies of pair creation, which is one of the 
possible purposes of the XFEL facilities, will require thor- 
ough knowledge of these characteristics and their depen- 
dence on the parameters of the laser field. 

In the present work, we address the problem of the 
pair creation in the tunneling regime. We propose a new 
variant of the Keldysh approach, the so called vicinal 
aprroximation, which takes into account important phys- 
ical properties of the system and simultaneously greatly 
simplifies the calculations required. We then find a com- 
plete analytical solution to this problem, deriving new 
analytical expressions for the angular distributions and 
energy spectra of the created leptons, as well as photon 
absorption spectra. The linear and circular laser polar- 
izations are studied in detail, and a brief outline of the 
general elliptic case is presented. Our preliminary results 
were discussed in Ref. |21i] . 



II. THEORETICAL FRAMEWORK 



Laser Parameters 
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The laser frequency cu is assumed to be small, so that 
it satisfies the adiabatic condition 



(2.1) 
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Unless stated otherwise we 
1, = a, e > 0, where — e 



where m the electron mass, 
use relativistic units H = c = 
is the electron charge. Secondly, the laser electric field 
strength £ is considered to be small in comparison with 
the QED critical field £c = m'^/e = 1.3 x lO^^ Ym-\ so 
that their ratio 



T = £/£c < 1. 



(2.2) 



Thirdly, the adiabatic (intensity) parameter of the prob- 
lem ^ is presumed to be large: 



e£ 



> 1. 



(2.3) 



Eas. (|2.ip - (|2.3p together describe the tunneling regime for 
the pair production. This regime may be achievable in 
future experiments, since it is expected |2,l22] that XFEL 
facilities will be able to produce electric field strengths 
up to £ ~ 10~^£c at frequency uj ~ 0.002m — 0.02r7i. 
Conditions (|2.ip - (|2.3p may be equivalently expressed by 
the combined inequality 



B. Amplitude and Probability 



(2.4) 



Consider the field produced by a plane electromagnetic 
(EM) wave and the Coulomb field of a bare nucleus in 
the nuclear rest frame, with the nucleus placed at the 
coordinate system origin. The nuclear Coulomb potential 
for an electron at radius r is thus —Za/r, where Z is the 
atomic number of the nucleus. 

The amplitude of a multiphoton process is given by a 
matrix element averaged over the period T = 2-k /uj of 
oscillations of the EM plane wave, see e.g. Ref. [13, [3] ■ 
Specifically for the case of the pair creation the amplitude 
reads 



{i:p_{t)\~Za/T\rr,M dt ■ (2.5) 



Here ilip_ and ijjp^ are the wavefunctions of the electron 
and positron. The subscripts ' — ' and '+' henceforth de- 
note electron and positron variables respectively, and the 
superscript * denotes a complex conjugate. Note that 
the electron and positron states are indexed by their mo- 
menta p± . The Floquet theorem guarantees that we may 
always write the wavefunction il>p^ of these particles in 
the form 



Both the electron and positron states are described here 
by complex-conjugated wave functions since they both 
represent the created particles. Together with the mo- 
menta p±, each wave function here is characterized also 
by the spinor index A± , which is not depicted explicitly. 

The probability of the multiphoton process is propor- 
tional to the square of the amplitude, |Mp^p _p which 
should be multiplied by appropriate, conventional statis- 
tical factors, see e.g. [13, [3 for details. Using this rule 
one writes the rate of the pair creation, W , as 



n— 1 



2^ \^^h+P- I (27r)3 (2^)3 



(2.8) 



Here n is the number of the photons absorbed, summa- 
tion over s± takes into account two spin states of each 
lepton s+ — ±1/2, while summation over A± refers to 
the spinor indexes (for the sake of simplicity this summa- 
tion is not presented in the following formulas explicitly). 
The Dirac delta-function imposes the energy conserva- 
tion law 



(2.9) 



In the case u) m, Eq. (|2.8p may be simplified by not- 
ing that the sum over n can be replaced by a correspond- 
ing integral, which gives the creation rate the asymptotic 
form 



W 



27r 



(27r)3 (27r)3 



(2.10) 



Taking partial summations and integrations here one 
finds expressions for the differential probabilities dW, 
which describe the energy spectra and angular distri- 
butions for the pair, as well as the photon absorption 
spectra. 



Volkov Wave Functions 



Presuming that the nuclear charge is not large, i.e. 



Za < 1 



(2.11) 



we treat the interaction of leptons with the Coulomb field 



V'p± (^, t) = 'fp± (*", t) exp {~iE±t) , (2.6) perturbatively to first order. The lepton state ijjp^ in this 



where E± is the mean quasienergy for the particle (see 
Section FlI Pp . and (pp^(r,t) is a periodic function of t 
with the period T. 

The matrix element (| |) in the integrand in Eq. (|2.5p 
means a conventional integral over spatial coordinates, 
that is 

(Vp_(oiiAiV';+(i)) = / rp.{r,t)rp4r,t)^ . (2.7) 



approach, which stems from Refs. 3, 12], can be described 
by a Volkov wavefunction that takes into account the 
lepton interaction only with the EM plane wave. 

Let A — A{^) be the 4-vector potential of the electro- 
magnetic plane wave, where the scalar $ = (kx), k being 
here the light 4-wavevector, fc^ = 0, and x = {t,x,y,z) 
is the usual spacetime coordinate. We use the conven- 
tional notation that = (a^), and a'^a^ = a^. In the 
Lorentz gauge, {dA) = 0, the state of the electron in the 
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EM plane wave is described by the Volkov wavefunction 
^_ (see e.g. Ref. [1^ for a derivation): 



1 - ^-^ijk){jA) ] exp (iS-) 



2 (kp) 



2e 



(2.12) 



where e > 0; p = (e,p) is the electron 4-monientuni (for 
clarity we temporarily drop the '±' subscripts for it), 
P = {Px,Py,Pz)', 7 = (7", 7) are the Dirac matrices; Up 
is a free Dirac 4-spinor, which satisfies 



ilP) Up = mup ; 



(2.13) 



and S- is the classical action for the electron in the elec- 
tromagnetic wave 



-{px) 



Pj^^{ipA) + lA^)d^. (2.14) 



Henceforth, we assume that the electromagnetic wave 
propagates along the z direction. Then k = (lj, 0, 0, w), 
so that (kp) = uj{e—pz), = Lu{t—z), and it is convenient 
to define 



1] ^ e -pz 



(2.15) 



Equation (|2.13p implies that p^ = m^, so that the leptons 
are created on the mass shell. It follows that 77 > and 

e= , (2.16) 



2-q 



2-q 



where = + • The 4- vector potential of the elec- 
tromagnetic wave may be written as 



A($) = ^a($) = ^ (0,a($; 



(2.17) 



where a($) is the polarization of the wave and we have 
chosen = 0. This polarization satisfies a($) • fc = 0. 
We may now write the electron wavefunction (|2.12p and 
action (I2.14|) in a clearer form: 



- (1 - ^(7'' - 7^)(7a)) e-P(*^-)^ ' (2- 



18) 



5_ = -{px) + ^ y"^' ^^[{pa) + jd* . (2.19) 



Similarly, the positron wave function and the action read 
*+ = (1 + ^(7° " l'){ia)^ °"P(^^+)^ ' (2-20) 
S+ = -ipx) ^ f^'^^^ [{pa) - "^a^y^ .(2.21) 



The spinor Up satisfies the Dirac equation (|2.13p describ- 
ing the propagation of a free positron, but in calculations 
it is convenient to express it conventionlly, via the spinor 
that describes an electron in the " lower" continuum using 
an operator of the charge conjugation Ctp = j'^ip*, 



up{ep) = -/^u*_p{-ep) 



(2.22) 



Here we have to specify an energy ±ep = ±(p^ -|- m^)^/^ 
of the state described by the corresponding spinor; the 
spinor 7ilp(— gp) describes an electron state with nega- 
tive energy, i. e. the state in the lower continuum. (For 
clarity spinors are presented below without reference to 
the energy of the state they describe.) 



D. Vicinal approximation 

The interaction of leptons with the Coulomb field must 
be strong, for otherwise the probability of the pair cre- 
ation would be very small. We therefore presume that 
the most important events during the pair creation take 
place in the close vicinity of the nucleus, i.e. we presume 
that the pair is created at distances which are compara- 
ble or smaller than the Compton radius rc = l/m. We 
will verify later that this presumption is correct by direct 
calculation (see Eqs. 13.261 and I4.26P . 

Having this presumption in mind we can presume now 
that \z\ < l/m in Eqs. ((TT5)) and ([^T^ . Importantly, 
in the adiabatic regime cj ^ m, so it follows that 



uj\z\ < 1. 



(2.23) 



We call inequality (j2.23p the vicinal approximation, since 
it arises from the close vicinity of pair creation to the 
nucleus. 

Since loz is small, we can expand the action (|2.19[) in 
a Taylor series about $ = to first order in loz, and 
replace argument uj{t — z) with ut in the pre-exponential 
term of the wave function (|2.18p . Thus in the vicinal 
approximation the Volkov wavefunction of the electron 
is greatly simplified: 



■i>{r,t)^Q{iut)^exp 



i[p{t)-r- j E{s)di 



where Q{ujt) is the matrix 



Q{u;t) = 1 



f^(7° 
2ry 



7')(7a(wt)) 



(2.24) 



(2.25) 



p(t) — {pxiPyiPL{t)^ is the quasimomentum and E{s) 
is the quasienergy. We call pl {t) the longitudinal quasi- 
momentum. Further applying Eqs. (j2.16p to the Taylor 
expansion permits us to present the longitudinal quasi- 
momentum and quasienergy as 



PL±{t) = 



E±{t) = 



-I- (k =F ^ma{u;t) ) — 
2^ 

+ (k =F ^ma{ujt)Y + 
2ry 



(2.26) 
(2.27) 



where n — {px,Py,0). Again for the sake of clarity we 
have suppressed the labels ± in all variables on the right- 
hand sides of these equations, except in front of ^ where 
it is essential. 
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The mean quasienergy E can then be derived from 
Ea. (|2.27p . via its definition 



E± = 



E±{s) ds. 



(2.28) 



Note that from Eq. (|2.24p we then may write ^{r,t) = 
e~'-^'(p(r, i) just as we did in Eq. (|2.6|) . where (p{r,t) is a 
periodic function of t with period T. 

The advantage of the vicinal approximation is that the 
time and space variables are decoupled in the wave func- 
tion (|2.24p . greatly simplifying evaluation of the integrals 
in the amplitude (|2.5|) . That is, the amplitude becomes 

Mp^p_ = ^ / IC{t) V{t) exp {iS{t)) dt, (2.29) 
^ Jo 

where the factor IC{t) includes the spinor variables 

(2.30) 

V{t) is the Fourier transform of the Coulomb potential 



V(t) 



AnZa 



(2.31) 



and S{t) is the contribution of the total quasienergy 
E{t) = E+{t) + E_{t) to the action 



S{t) = J E{s)ds 



(2.32) 



in which 



A = 



f V- 
1 K^rj^ — K^rj^ 



1 



V- 



2 f^+V- 

X = — I m + H 



- (mm'- 



(3.3) 
(3.4) 

1/2 
I ■ 

(3.5) 



From Eq. (|3.2p . the mean total quasienergy is clearly 



E = X{l/2 + p' + x'). 



(3.6) 



The mean quasienergy has a minimum when Kj^ = /t_ = 
0, rj+ ^ rj^ — m(l + ^ 12)^1"^ . At this minimum 



E = 2m(l + i^l^fl'^ ~ V2mC 



(3.7) 



in the tunneling regime. This minimum is the threshold 
energy required for pair creation, and it grows linearly 

with \: E = 2m(l + ^2/2)^/2 ~ ^2™^. 

Also from Eq. (|3.2p one finds the action (|2.32p is 



S{t) 



m 



EX, 

— uot — - — cosujt (sinwt + 4p) 

TO 2m 



(3.8) 



Eq. (|2.ip ensures that the coefficient m/uj in front of 
this expression is large. Consequently, the integral over 
time in Eq. (|2.29p can be evaluated by the saddle point 
method. Labeling the saddle points as tj, j = 1,2, . . ., 
we find for the amplitude 



p+p- 



27r 



iE'{t,) 



1/2 



(3.9) 



III. LINEAR POLARIZATION 

Consider an electromagnetic plane wave with linear po- 
larization 



a($) = (0,- sin $,0,0). 



(3.1) 



In this section we write the lepton momenta as p± = 
(e±, K^, K^jPzi) = (£±,p±). The £ and H. superscripts 
denote the electric and magnetic field directions respec- 
tively, which correspond to the x and y coordinate direc- 
tions. Further, k± = is the lepton momentum 
transverse to the wavevector. 



A. Probability 

For the polarization given by Eq. (|3.ip the total 
quasienergy found from Eq. (|2.27p may be rewritten in 
the simple form 

E{t) = A ( (sin ojt + + x^) , (3.2) 



where the prime denotes differentiation with respect to 
t. 

The saddle points must satisfy the following condition: 



S'(t^) = Eit,) = 



(3.10) 



In other words, the pair is produced at the moment it 
has zero quasienergy. From Eq. (|3.2p this condition is 
equivalent to 



sin ujt 



-(3±ix. 



(3.11) 



We are interested only in the saddle points with real part 
of tj lying in [0, T] , so there are only two pairs of complex 
solutions of Eq. ()3.1ip : one solution of each pair lies in 
the upper-half complex plane and another in the lower- 
half complex plane. Since pair creation requires absorp- 
tion of ener gy f rom the laser field, the theory of adiabatic 
transitions [23 . [25l | specifies that only the saddle points 
lying in the upper-half complex plane contribute to pair 
creation. Hence we have two saddle points 



w^i = sin ^(/3 + zx) + 7r 



Ujt2 



sm 



\-P + ^x) 



mod 2iT 
mod 27r 



(3.12) 
(3.13) 
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These saddle points have the properties that 

= 3(wt2) > (3.14) 
5R(wt2) =7r-5R(wti) mod27r. (3.15) 

Examination of Eq. (|3.8|) reveals that the action S{tj) 
consequently satisfies 



5[5(<i)] -3[5(t2)] >0 
Ett 
to 



n [S{t2 



^[S{ti)] mod27r. 



(3.16) 
(3.17) 



The pre-exponential terms of the sum in Eq. (|3.9[) 
prove be same for either saddle point. This fact combined 
with Eqs. ()3.16p and ()3.17p permits us to express the 
amplitude in terms of only one saddle point, for example 
h. By substituting Eqs. ((3T6)) and ((3l7| into Eq. 
the probability becomes 



TT \E'{h)\ 






1 1 + cos 









exp(-23[5(ii) 



(3.18) 

The phase within the cosine is large, being cx m/w, so 
that as the lepton momenta vary slightly we expect the 
cosine term to oscillate very rapidly. As a result, we 
may neglect the contribution of this term to the integral 
in Eq. (|2.10p . and hence to any of the spectra or rates 
discussed in the present work. We therefore ignore this 
term henceforth. 



B. Pair Creation Rate 



It now follows from Eq. p.lSp that the pair creation 
rate (PTTU)) is 

'"^j^J \E'{t,)\ (27r)3(27r)3- 

(3.19) 

Let us find the minimum of 3[iS(ti)] as a function of 
electron and positron momenta. For this purpose, it is 
convenient to define P = (k^, k^, r/^, k^, k^, rj-) and to 
write all functions as explicit functions of both P and 
noting that the saddle point ti ~ ti{P). The minimum 
of3[5(ti,P)] then satisfies Vp3[5(to, -Fb)] = and cor- 
responds to the most probable configuration of electron 
and positron momenta. 

It is instructive to compare the current problem to the 
well-known static case for homogeneous electric and mag- 
netic fields, where the fields are orthogonal and equal 
in magnitude: i.e. E ■ B = 0, = E^. In this 
case any charged particle would accelerate mostly in the 
E X B direction regardless of the sign of its charge. We 
may expect therefore that the minimum of 5[iS(ii, P)] is 



achieved when the electron and positron move along the 
z-direction. 

By considering the partial derivatives of 3[5(io,-P)] 
with respect to and k^, it is straightforward to verify 
that at the minimum of 3[5(ti, P)], 



= K = 0. 



(3.20) 



Further, from the symmetry of the problem, we expect 
that the electron and positron momenta should be equal 
at the minimum of 3[5(ii,P)], so we expect 77+ = 77-, 
since rj defines the z-direction momentum pz by Eqs. 
(I2.16p . This condition together with Eq. (|3.20p permits 
the form of the quasienergy (|2.27p to be simplified, and 
so the action p.32p becomes 

S{tuv) = ljj' {"i'[l + C'sin2(cji')] dt', (3.21) 

at the minimum of 3[5(ti, P)]. 

Applying Eq. (jX^ to Eqs. (EH), (ES]), and (IXTTj) . we 
also find that /3 = 0, x = {w? -\-rfY/'^ /m£, and s\n{ujti) = 
ix respectively. We presume that at the minimum rj ^ m 
(we verify this in Eq. (|3.24p below). Then for ^ ^ 1, we 
find 



1 



M = - 1 + ^ 



2 \ 1/2 



< 1. 



(3.22) 



It follows from the inequality that sin{ujt') ~ ujf in Eq. 
([3^ . After applying Eq. ((3?22)) to the action ((3?2T|) we 
then have 



2 (m2+ 77^)3/2 
3 J- m?!] 



(3.23) 



The imaginary part of this expression is minimal at 77 
so Q[5(ti,P)] is minimal at 



„ , 771 777 

P„H 0,0,-^,0,0,-^ 



with minimum 



^[S[h,Po)]^^. 



(3.24) 



(3.25) 



It is important to note that from Eq. (|2.26p the quasi- 
momentum at the minimum is 



(3.26) 



Equation (j3.26p shows that the typical momenta in the 
Coulomb matrix element of Eq. (|2.5p are |p| ^ m. This 
indicates that the distance from the nucleus within which 
the pair is created is comparable to the Compton radius. 
I.e. r < l/\p\ ~ 1/771. This estimate supports the validity 
of the vicinal approximation introduced in Section III Dl 
We may expand 3[iS(ti, P)] in a Taylor series in powers 
of P about the minimum Pq. This expansion may be 
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calculated by using Eqs. ([S^-jSll) and dSU). The 
result to second order in P is 



(3.27) 



where Hp^ is the 6x6 Hessian matrix of second order 
partial derivatives of 5[5(to, P)] a± P — Pq and {P—Pq)^ 
denotes the transpose of (P — Po)- In the regime ^ ^ 1, 
this matrix is 



Hp. 



V3 
2m?T 



1 















0' 







2 






















7 

3 








1 

3 


1 


3 

w 






























2 













1 

3 








7 

3- 



(3.28) 



The higher order 1/^^ terms are included in the matrix 
so Hpg is not singular. From Eqs. p.27p and (|3.28p we 
may thus write for ^ ^ 1 



A{P) = 



12to2 



2m^ 



£ \2 



+ 



(Ki + (a«^ - 4) 



Am? 



8^2 m2 



(3.29) 



where (5?7± = ??± — m/V2 is the variation of ?7±. 

The factor v^/J^ in Eq. ((3:27| is large since <C 1, 
so we can apply the saddle point method in order to 
evaluate the integral in the pair creation rate (|3.19p . We 
now change integration variables in the integral p.l9p 
from <f'p-(Pp+ to d^P, so the minimum point Pq is the 
only saddle point for the new integral. The differential 
pair creation rate thus becomes 



where the coefficient 



\E'{ti)\ \di] 



(2^) 



exp 



(3.30) 



2x/3 



(3.31) 

and the saddle point ti is evaluated at Pq. 

The only spin dependent term in B is IC{ti). In Ap- 
pendix 1^ we verify that the spin dependent factors give 
the following contribution to the creation rate 
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(3.32) 



Further, Eqs. p.26p and (|2.3ip permit us to calculate 
the matrix element of the Coloumb potential 



\Vih)\' 



2nZa 



(3.33) 



and Eqs. ([3?22)) . and ((3?24)) together produce 

^'(^i ) = i2V5J^m'^ in the ^ > 1 regime. From Eq. 
(I2.16p . the partial derivative dpz/drj — —3/2 at the sad- 
dle point. Hence, 



B 



8^2 z^a^u; 



exp 



(3.34) 



V3 

Straightforward integration of Eq. p.30p produces the 
following result for the pair creation rate in the linear 
polarization case: 

3 



W 



7^^/2 V2%/3 



exp 



2^/3 



(3.35) 



We will check later (see Eq. 14.371) that this formula agrees 
with the previously known results. 



C. Positron Spectrum 

Let us find the energy spectrum of the positrons. We 
may integrate the differential pair creation rate in Eq. 
(|3.30p over the electron momenta dn^dK^dr]^, so that 



dW = W+dn%dHi'^d7i 



(3.36) 



Clearly W+ represents the differential positron creation 
rate, as it is only a function of the positron momenta 
(k^, K^, 77-i_) . We then have 



V3 

W+ = B+ exp ( - —A+{K%,K^,fi^ 



(3.37) 



where Q = and the exponent and coefficient are 

respectively 



(3.38) 



4 Z^a^ f T i 2^\ 

+ ~ (2^)5/2 V7 UVsj ^""^y T )■ 

(3.39) 

The positron energy spectrum is dWjdE^, where i?+ is 
the positron mean quasienergy in the laser field. (An 
empirically important property of the created positrons 
is their energy outside the laser field. If one assumes that 
the positron leaves the focus of the field adiabatically, 
then its energy outside the laser field is exactly its mean 
quasienergy i?+.) 

In order to find the positron energy spectrum, we must 
write dW (|3.36p as a differential form which includes 
so that a change of variables is required. The trans- 
formation of dW under a change of variables vijj) = x 
from X = (xi, . . . , Xn) to y = (j/i, . . . , y„) for the saddle 
point = i'(yo) reads 

dW = B+\jy.v\ exp({y~yQf{Jy„i'fH+{Jy„i'){y-yo)jd''y 

(3.40) 
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where JyaV is the Jacobian matrix of v evahiated at sad- 
dle point ?/oj and is its determinant. From Eq. 
p.38p the matrix iJ+ is simply the diagonal 3x3 matrix 



V3_ 



diag{l,^^8f/7}. (3.41) 



From Eqs. jMZl), dMS]), and ([33]) and the fact that 
77 > 0, one finds that 



v+^E+~ [El -K^- m\i + em] 



1/2 



(3.42) 



Equation ()3.42p defines the required coordinate transfor- 
mation ve{h%th'^,E+) = («+,«+,??+)• At the saddle 
point, E+o = to(3 + ^2)/2V2 ~ 1'^^, f or g > 1. 
Substituting Ve and the saddle point into Eq. (|3.40p and 
then integrating over and k^, one finds the positron 
energy spectrum 



dW 



(. 



dE+ (2^)3/2\/7 e V2\/3 

32V3 f E+ 
'7^ \ m " 



5/2 



exp 



2V3 



X exp 



2V2 



(3.43) 



D. Angular Distribution 

Let us find the positron angular distribution, dW/ dQ^ 
which indicates the direction of the mean velocity of the 
positron within the focus of the laser field. (If param- 
eters of the field are known, one can extract from this 
information the angular distributions of the leptons out- 
side the laser field.) We employ the spherical coordinates 
(Ip+1,6', (p) and dfJ — sm9d9d(p. In this coordinate sys- 
tem, 



|p+| cos ip sin 6 , 
\p+ \ sin (/3 sin , 



(3.44) 



V+ = {m' + \p+fy^'' -\p+\cose 



which follows from Eq. (|2.16p . These together define 
the transformation i/Q{\p^\,9,ip) — (k' 



The 



saddle point in the spherical coordinates is located at 
= and 6 = 0. Note that the azimuthal 

angle ip is not well-defined at this saddle point, so we 
keep all orders of ip in the differential positron creation 
rate exponent A+. After application of Eq. p.40p and 
integration over one finds 



dW 



1 Z^a^m 



(—)' 

dn 167r2V2 C \2V3J 



exp 



2x/3 



X exp(^-^{cos'^ p + ^siri^ p)0^y (3.45) 

We assume in this derivation that the parameter Q = 

is small, i.e. t/ ^ 1, so that the angular distribution is 



integrable. The condition Q 1 represents a special 
case of the tunneling regime, and it is equivalent to the 
additional inequality m/uj S> ^ (see Eq. ()2.4p ). The 
operational parameters of an XFEL may fall just inside 
this regime 



E. Photon Absorption Spectrum 

Let us find the spectrum of the number of absorbed 
photons in the pair creation process. In Section III Bl we 
used the adiabatic condition uj <^ m (|2.ip to remove the 
sum over n from the pair creation rate in Eq. (|2.8p . We 
now define the photon absorption spectrum as 



|AWP<5(i?-n.)|gt|^, (3.46) 



where the total mean quasienergy E = E^ + 
Note that by comparison to the pair creation rate (|2.8p . 
Y,^ Wr, = ly, or / Wndn = for w < to. 

We now replace the Dirac delta function by its Fourier 
transform: 



6{E — nuj) = — / e" 
27r 



{E-r 



'dr. 



(3.47) 



From Eq. p.l8p we then write 

l^(ii)V(ii)P 



E 



\E'{h)\ 



e.p{^rE-2^[Sit,)]}^^dr, (3.48) 



Just as shown in Section IIIIBI for the integral in Eq. 
(I3.19p . the integral over momenta in Eq. (|3.48p may be 
evaluated by the saddle point method. From Eqs. (|3.2p - 
p.5p . the irE term in the exponent of Eq. p.48p does 
not have a 1/^ factor, so that when ^ 1 the saddle 
point Po is the same as it was in Eq. p.24p . Upon 
application of the saddle point method, the irE factor is 
instead expanded about Pq to first order. Hence 



27r 



where the exponent 



exp[C(P,r)]-— gdr. 



(2^)6 



(3.49) 



CiP,T) = irV pE{Po) ■ {P ~ Po) 



2V3 



A{P). (3.50) 



With reference to the definition of A p.27p and employ- 
ing some linear algebra, it can be shown that 



W„ = W— I exp ( iT{EiPo) - nu;) 



exp( -^[\/pE{Po)fHpypE{Po))dT. 



(3.51) 



8 



From Eqs. ([331) -(ES]), (jSTM)) . and ((3:28)) we then have 
that for ^ > 1 



Wn = W— / exp 



V3, 



me 

V2 



X exp 



J^m^^'^T^ 1 dr. 



128 



(3.52) 



Evaluating this integral produces the photon absorption 
spectrum in the tunneling regime: 



5/2 



(2^)3/2 V3 e Uvsj """H y 



X exp 



2 \ 2 



;V2 



(3.53) 



IV. CIRCULAR POLARIZATION 

Consider an electromagnetic plane wave with circular 
polarization 



a(<f>) = (O, - sin $, cos $, O) . 



(4.1) 



In this section it is natural to employ the cylindrical 
coordinate system, so we write the lepton momenta as 
p — (e, Kcos^s, Ksin(^,p^) — (e,p), and further k = 
(k cos iy9, K sin (p). Here the subscripts ± have been omit- 
ted for clarity. It is also convenient to write the mean 
lepton angle and difference of lepton angles respectively 
as 



ip = and = 



A. Probability 



(4.2) 



For the polarization given by Eq. (|4.ip . the total 
quasienergy found from Eq. (|2.27p may be written as 

E{t) = a(x + /3sin(c^i - - a)) , (4.3) 

where A is the same as in Eq. p.3p . and here 
2 |77_/t+ - ?7+K-| 



/3 = 



1 2/1 , ci\ , , K^??- + «:-»7+ 



m2^2 



(7 = tan I tan < 



(4.4) 

(4.5) 
(4.6) 



From Eq. (|4.3p the mean total quasienergy of the lep- 
ton pair is 



(4.7) 



The mean total quasienergy is minimal at k+ = k_ = 
and 77+ = 77- = m(l -I- f^)^/^, with minimum value 



E = 2m(l +^2)1/2 _ 



(4.8) 



in the tunneling regime. This minimum is the threshold 
energy required for pair creation, and just as in the linear 
polarisation case, it grows linearly with ^. 

It also follows from Eq. (|4.3p that the action (|2.32p is 



„, . m f E A/3 , , 

6(<) = — — ujt H cos(cjt — (/3 — fj) 

Lu \ m m 



(4.9) 



Just as in Section IlIII the action has a coefhcient m/uj, 
which is large by Eq. (|2.ip . Hence the integral over 
time in Eq. (|2.29p may be evaluated by the saddle point 
method. For saddle points tj the amplitude Mp^p_ has 
the same form as in Eq. (|3.9p . 

The saddle points must satisfy E[tj) — S'{tj) — 0, 
which here is equivalent to 



X 

s\n{Lotj — Lp — a) = — — . 

P 



(4.10) 



The right side of Eq. (|4.10p is real and we presume 
|x//3| > 1 (verified in Eq. 14.201 below), so there is only a 
single pair of saddle points with real part in [0,r]: one 
lies in the upper half complex plane and the other in the 
lower half complex plane. Just as in the linear polariza- 
tion case we choose the saddle point lying in the upper 
half complex plane [13, HI], and label it to. With refer- 
ence to Eq. (|3.9p . since there is only a single saddle point 
the probability becomes 



|M„ 



|/C(to)pV(to)^ 



27r 



\E'{h)\ 



exp I 



23 [S{h)\ 



(4.11) 



Pair Creation Rate 



From Eq. (|4.1ip the pair creation rate (|2.10p is 



W 



\lE'{t^)\ 



(27r)3 (27r) 



3 " 



(4.12) 

This formula differs by a factor of 2 from the rate in the 
linear polarization case (|3.19p since there is only a single 
saddle point here. 

Let us find the minimum of 3[iS(to)] as a function of 
the electron and positron momenta. For this purpose, 
it is convenient to define P = 77- , 0, (/j) and 

to write all functions as explicit functions of both P and 
t, noting that the time saddle point to = to{P)- The 
minimum of 3[5(to,P)] satifies Vp9[5(to, ^')] = 0. 

The circularly polarized EM plane wave imparts angu- 
lar momentum to the lepton pair, so we expect that k+ 
and K_ are not zero at the minimum. By symmetry, we 
expect that at the minimum the electron and positron 
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momenta are the same in size and the total momentum 
of the pair is conserved. Hence we deduce that at the 
minimum: 



V+ = V- = V- 



(4.13) 
(4.14) 



Eq. (|4.13p implies = k_ = k. Moreover, if we 
choose the azimuthal angles so that (p_ € [0, 27t] and 
tp+ € [tTjStt], then this condition means that 



(4.15) 



It follows by Eqs. ()4.2p that at a minimum Po of 
3[5(io,P)] ■ 



I , and 



TT 57r 
2' T 



(4.16) 
(4.17) 



Applying conditions to Eqs. ((i3)) . and 

at this minimum, we also find that /? = 2K/m^, x — 

[m2(l + + ,y2 _|_ ^2j^^2^2 ^^^^j ^ ^ ±7r/2. The ± sign 

here is due to the fact that the angle cr is not well-defined 
under conditions (|4.13|) . The total quasienergy (|4.3p is 
then simplified, so that the action (|2.32p becomes 

1 /■*" 

5(to) = - / [m'^{l+S,'^)+r]'^+K'^±2mK^cos{ujt'~ip)]dt' 

^ . (4-18) 
at the minimum of 5[5(io, ^')]- 

Consider now the partial derivatives with respect to 
rj and k of 5[5(io)] (|4.18p . which must be zero at the 
minimum. This produces two further equations, which 
may be divided by one another to derive the relation 



m2(l + f )-7/2 



(4.19) 



Just as in the linear polarization case (Section IIII Bp 
we presume ~ m. In fact, since the circularly polar- 
ized wave is equivalent to two distinct linearly polarized 
waves with particular phase and polarisation, and since 
by its definition (|2.15p ur] is a frame invariant quan- 
tity, we expect ?/ = m/y/2 at the minimum. (This is 
verified below in Eq. 14.241 ) Then from Eq. (|4.19p . 
K = m(l/2 + ^2)1/2 ^ jn^^ for C > 1- 

At the minimum Eq. (|4.10p now becomes 

from which it can be seen clearly that |x//3| > 1. For 
^ » 1 we find 



Mo-^l = ^(i + ^j «1 



(4.21) 



It follows from the inequality in Eq. (|4.2ip that 

cos(wi' - 0) ~ T [1 - i^^t' - 0) V2] , (4.22) 



in the integral in Eq. (|4.18p . Applying Eq. (I4.2ip to the 
action in Eq. (|4.18p we find that 



2 (m2+ ,72)3/2 ^ 

^ / ^ , W(l+2e2)+?72 



mP"q 



CUT] ' 



(4.23) 



The imaginary part of this expression (Eq. (|3.23p ) is min- 
imal at 77 = m/^/2. Importantly, this minimum is in- 
dependent of the mean angle ip, so the minimum Pq is 
degenerate in ip. This can be verified by noting the par- 
tial derivative of the action (|4.9p with respect to (p is 
purely real. Hence, !3[5(to,^')] is minimal at 



Pq 



m m TT 



2'T 



(4.24) 



with minimum value 



3[5(to,Po)] = ^ 
The quasimomentum (j2.26p at the minimum is 



cos Lp, sin Lp, - 



2V2J 



(4.25) 



(4.26) 



so in the circular polarisation the typical momenta in 
the Coulomb matrix element of Eq. (|2.5p is \p\ ^ m^. 
This implies that for ^ 3> 1, the distance from the nu- 
cleus within which the lepton pair is created is much less 
than the Compton radius. I.e. r < l/\p\ <^ 1/m. This 
estimate strongly supports the validity of the vicinal ap- 
proximation introduced in Section [II Dl 

From Eqs. (gS]) and (|¥TU)) . 5[5(io, P)] is independent 
of ip, so it is convenient to define P = {R, tp) and write 
the imaginary part of the action as 3[5(to,-R)]- We now 
seek to expand 9[5(to,i?)] in a Taylor series about Rq. 
This expansion is calculated using Eqs. ()4.3p - ()4.9p and 
(|4.10p and noting 5 [to] is a function of Rq. To second 
order in R abut Rq, the result of this expansion is 



^S{to,R)] = ^ + -{R-RofHR^{R-Ro) 



^[l+AR)], 



(4.27) 



where is the 5x5 Hessian matrix of second order 
partial derivatives of 3 [5(^0, R)] evaluated at Rq. In the 
regime ^ > 1, 



H 



Ro 



V3 



"1 











■ 





7 

6 





1 

6 











1 











1 

6 





7 

6 





.0 











2m2e2 



(4.28) 



From Eqs. (jOT)) and (g^S]), we then have 
A{R) 



7{5if_ + 5f]\) + 26ri^5rj+ ^ Sk\ + Sk^_ , ^2 r 2 



12m2 



2m? 



(4.29) 
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where Sr] = r] — m/-\/2, 5k = k — m^, and 5(j) = 4> ~ 7r/2. 

The factor in Eq. (|07|l is large, so we may 

immediately apply the saddle point method to evaluate 
the integral in the pair creation rate (|4.12p . Clearly the 
saddle point is the minimum Pq. The degeneracy of the 
saddle point Pq in its ip component means that its con- 
tributions includes an integral over ip from tt to 57r/2. 
However, all terms in the integrand of Eq. (|4.12p prove 
to be independent of tp, as expected by the cylindrical 
symmetry of the problem, so this integral merely contr- 
butes a factor of 27r. Thus we have the differential pair 
creation rate 



(4.30) 



where the coefficient 



\JC{to)V{to) \^ (dp, 
\E'{h)\ 



[dp.' 


2 

) 










Ro 


d{(t),ip) 



_ 2y3 

e ^ 



(4.31) 

and the saddle point to is evaluated at Rq. 

The only spin dependent term in B is JC{to). In Ap- 
pendix 1^ we verify that 



(4.32) 



will then be the average over one period of the plane 
wave. We thus have linear pair creation rate 



W 



2V3 



5/2 



X exp 



T\ cos(a;t)| 



dt. 



(4.36) 



This can be evaluated by the saddle point method, since 
<C 1, with saddle points at t = 0, T . Hence 



W 



2 Z'^o?m, 


( T 


T 2^ 


( 




^V2 





5/2 ^oo 



exp 



cxp 



2jf 



~y 2~ 



dt 



(4.37) 



which is precisely the result derived above (|3.35p . The 
linear polarization pair creation rate p.35p is therefore 
also in agreement with the static crossed field or homo- 
geneous field results. 

The form of the circular polarization pair creation rate 
(I4.35P agrees with the form of the result derived without 
use of the vicinal approximation [9], and both rates are 
also in exact agreement with the result derived by use of 
the polarization operator [23| . 



From Eqs. ([23T|l . ijiTS]) . and (jiTM]) one finds the 
Coulomb potential is 



C. Positron Spectrum and Angular Distribution 



|V(io)|^ 



2'KZa 



(4.33) 



Equations (|4.3p and (|4.2ip in the regime C ^ 1 produce 
E'{tn) = ^^2^/3J^m^ and by Eq. ((2T6l) . dp^Jdr^ = -(^ at 
i?o. Further, the determinant \d[ip+,ip^) / d{(p,ip)\ — 2. 
We then have, see Eq.dHSlI), 



„ 8^2 ^2^2^ 

B = — 1= — = — ^ exp 



(4.34) 



Straightforward integration of Eq. (|4.30p produces the 
following result for the pair creation rate in the circular 
polarization case within the tunneling regime: 



W 



20r V2%/3 



5/2 



exp 



2V3 



(4.35) 



This result is precisely the pair creation rate for a static 
crossed field combined with a Coulomb field [2^, 
which is expected as the electric field strength is con- 
stant in a circularly polarized plane wave. The result is 
also in agreement with the pair creation rate for a static 
homogeneous field combined with a Coulomb field [28j . 
The pair creation rate for a linearly polarized EM plane 
wave may be derived from this result by assuming that 
the electric field in Eq. (|4.35p can be considered as time- 
dependent, E = £co8{ujt). The linear pair creation rate 



Let us find the spectra of the positrons. We employ the 
same approach to find these spectra as done for the linear 
polarization case in Section llll CI In order to find the dif- 
ferential positron creation rate W+ we integrate the dif- 
ferential pair creation rate in Eq. (|4.30p over dK-dr]-d<p, 
and write 



dW = W+dK+d'n+. 



(4.38) 



Note that by definition (|4.2p the angle (p describes the rel- 
ative orientation of the electron and positron momenta, 
so we must integrate over this variable. We then have 



a/3 

W+=B+e^p[ - —A+iK+,ri+) 



(4.39) 



where the exponent and coefficient are respectively 



(2^)5/2^7 m \2Vi) 



B. 



exp 



(4.40) 
(4.41) 



We may write down the transformation of dW under 
a general change of variable, as was done in Eq. (I3.40p . 
Such a change of variable is required to find the positron 
energy spectrum or the positron angular distribution. 
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From Eq. (|4.40p the matrix for the circular polariza- 
tion case is the diagonal 2x2 matrix 



x/3 



diag{ 1,8/7}. 



(4.42) 



The positron energy spectrum is dW/dE+. From Eqs. 
(I22I1), (P^ . and (P7T|1 and the fact that 77 > 0, we have 



E+- [El- K\~m^{l+^'^ 



,1/2 



(4.43) 



This defines the transformation ve{h+, E^) = [k,^,?]^), 
which permits dW to be expressed in terms of the dif- 
ferential dE+. At the saddle point Pi, E+q ~ m{3 + 
4^2)/2V2 ~ %/2m^2^ for ^ > 1. Substituting the co- 
ordinate transformation i'e and the saddle point into 
Eq. (|3.40p and then integrating over k+, one finds the 
positron energy spectrum 



1 



dW _ 

dE+ ^ 27r\/7 \.2V3j 
2V3 fE^ 



cxp 



X exp 



ev2 



2V3 
2-1 



(4.44) 



Due to the cylindrical symmetry in the circular po- 
larization case, the angular distribution here is simply 
dW/dO, where 9 is the polar angle. In the spherical co- 
ordinates (|p+|, ^) 



IP 



= IP+lsine* 



I cos 6 



(4.45) 



which follows from Eq. (|2.16p . These relations define 
the coordinate transform ve{\pj^\,6) = (k_|_, ry_|_), and the 
saddle point in spherical coordinates is located at | ~ 
TO^^/2, e = V2/C for C > 1. After application of Eq. 
(I3.40p and integration over |p+ 1 , we have 



dW Z'^c?'m^( T 



de (27r)\/7 V2\/3 
7T 



cxp 



X exp 



(4.46) 



D. Photon Absorption Spectrum 

Just as in the linear polarization case, the photon ab- 
sorption spectrum in the circular polarization case is (see 
Eq. ^M) 



Wn = W— I exp ( iTiE{Ro) ~ nu] 



exp (^-^[V rE{Ro)Y H^l ^nE{Ro)^dT. 

(4.47) 



From Eqs. ([SJl), gU-gJ]), (liTM)) . and (liTM)) one then 
finds that for ^ 3> 1 

Wn = J exp (^iT{2V2m^^ - nw)^ 



X exp 



-J^m^^V^ I dr. 



(4.48) 



Evaluating this integral produces the photon absorption 
spectrum: 



(4.49) 



X exp 



OJ 



V. ELLIPTIC POLARIZATION 

Consider an EM wave with the elliptic polarization 

a = (0, - sin$, 6 cos$) . (5.1) 

The parameter 6, which is presumed to satisfy < 6 < 1, 
measures the ellipticity: 6 = gives the linear polariza- 
tion Ea. (|3.ip . b — 1 corresponds to the circular polar- 
ization Eq. (|4.ip . Wc verified that the probability of the 
pair creation achieves its maximum when the electron 
and positron momenta satisfy 



V- = v+ = v 



(5.2) 
(5.3) 



We already know that this condition holds for the linear 
and circular polarizations, as is discussed in and after 
Ea. p.20p . and Eq. (|4.13p . Further analyses reveals that 



K = 6 + 1/2)^^^TO ~ b£,m 



V = 



V2 



(5.4) 
(5.5) 



Here k = |/«_| = Calculating the behavior of the 

imaginary part of the action Eq. ()2.32p considered a func- 
tion of all momenta in the vicinity of its minimum we 
find 



AiP) 



V3 



[1 + -4(P)] , 



(5.6) 



2tv? 



+ 



Am? 



8^2 



(5.7) 



Here deviations of momenta from the values, which sat- 
isfy Eqs. ((0)l .( (53|) is implied. 

Calculations, which are similar to the ones discussed 
previously for the linear and circular cases give the to- 
tal rate Wdi of the pair production for a general elliptic 
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polarization. It can be conveniently presented using a 
coefFecient k^u, which distinguishes it from a rate for the 
linear polarization Wun, which is given in Eq. p.35p 



lin 



■Coll 



(1-62)1/2 



erf 



1/2 



(5.8) 
(5.9) 



Here erf(a;) is a conventional error- function 

2 r 

eii{x) = ^= I exp(-z^)dz . (5.10) 

For h — Q one finds that /cgii — 1 (since T ^ 1), 
which means that Eq. (|5.8|) correctly reproduces the rate 
for the linear polarization. For 6 — + 1, Eq. (|5.9p gives 
fccii = (ttVS/J^)^/^. Substituting this result in Egs. fO)) . 
(|3.35p one reproduces the rate for the circular polariza- 
tion Eq. (l4.35|) . The elliptic polarization was discussed 
previously by Milstein et al. "2(1| presuming that the 
polarization is not close to the circular polarization, i.e. 
b deviates significantly from 1. In that case, the error- 
function is close to unity and fcdi ~ l/Vl — 6^, which 
agrees with Ref. 20] . An advantage of Eq. (|5.8p is that 



it describes any polarization, without any restrictions. 

Eqs. (|5.6[ ). (|5.7p allow one to consider the spectral and 
angular distributions of the lepton pair for the elliptic 
polarization, though we will not dwell on this issue in 
this work. 



VI. DISCUSSION 

A. Above-Threshold Pair Creation 

The positron energy and photon absorption spectra are 
in essence non-normalised probability distribution func- 
tions, which describe the probability of pair creation as 
a function of positron energy and number of absorbed 
photons respectively. The asymptotic expressions ob- 
tained for the positron energy spectra (Eqs. (|3.43p and 
(I4.44P ') and photon absorption spectra (Eqs. p.53p and 
(14.49^ 1 are all in Gaussian form. Since the maximum of a 
Gaussian probability distribution function is its expecta- 
tion value, we may immediately write down the expected 
positron energies and photon absorption numbers, which 
are 



{E+)lin 

{n)un 



2V2' 



{E+),,rc - V2me, (6.1) 
{n)circ = , (6.2) 



where the subscripts 'lin' and 'circ' denote the linear and 
circular polarization respectively. 

The threshold quasienergy Eq required for pair cre- 
ation was derived in Eqs. (|3.7p and (14. 8p for each po- 
larization case. By symmetry, i?o ~ so we have 
that 



E. 



+0i„ - 



71' 



E. 



+0c 



m^. (6.3) 



Similarly, the energy conservation law E = nui (see 
Eq. ()2.9p ) provides that the threshold photon number is 
riQ — 2Eq/uj. In the ^ ^ 1 regime, the ratios of the 
threshold and expected values are then 



(n) _ 

'^0 lin E+Q 



in) 



{E, 



no arc E+o ^^^^ 



2' 



(6.4) 
(6.5) 



The result for the circular polarization is in agreement 
with previous results derived without use of the vicinal 
approximation [5, 9], and we have now shown that a sim- 
ilar result holds for the linear polarization case. 

The positron energy spectra and photon absorption 
spectra are broad, due to the 1/!F^^ and uP' iTrv?^^^ fac- 
tors in their respective exponents. However, compared to 
the typical scales of the spectra, defined by the expected 
values and (n), the spectra are actually quite nar- 

row in the tunneling regime. That is, from Eqs. (j3.43p 
and (|4.44p . the width of the energy spectra is 



E. 



(E^ 



(6.6) 



and similarly from (|3.53p and (|4.49p . the width of the 
photon absorption spectra is 



6{n) - VT—C^ < — 



(6.7) 



By Eqs. (|6.4p and (|6.5p . in the ^ ^ 1 regime the thresh- 
old energy £'+0 ^ (-^+) and threshold number tt-q ^ (n). 
The narrow width of the spectra in comparison to the 
their typical scales means that the spectra are strongly 
suppressed at the threshold. In other words, only lep- 
tons (absorbed photons) of energies (number) a factor 
of ^ above the threshold value contribute significantly 
to the pair creation rate. By arialogy to the (tunnel- 
ing) multiphoton ionisation case |19|, this effect is called 
above-threshold pair creation. This is shown in Figs. [1] 
and[2l Note that in these plots we choose uj ~ IkeV, and 
T = 0.02, 0.03 based on possible operational parameters 
of XFELs 0. 



B. Angular Distributions 



The angular distribution (|3.45p in the linear polarisa- 
tion case is a Gaussian function of the polar angle 6, so 
that the expected direction of positron emission is along 
the 9 = direction. Moreover, the factor of the sin^ ip 
term in the exponent means that the the angular dis- 
tribution is strongly suppressed in the direction of the 
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FIG. 1: Positron energy spectra (normalized to unity at their 
maxima) for (a) the case of linear polarization, and (b) the 
case of circular polarization. Here ui = 0.002m with T = 0.02 
(solid line) or .7^ = 0.03 (dashed line), so that ^ = 10 or 
15 respectively. Arrow A (B) indicates the location of the 
threshold energy E+o for T = 0.02 (0.03). 



FIG. 2: Photon absorption spectra (normalized to unity at 
their maxima) for (a) the case of linear polarization, and (b) 
the case of circular polarization. Here u = 0.002m with JT = 
0.02 (solid line) or ^ = 0.03 (dashed line), so that ^ = 10 
or 15 respectively. Arrow A (B) indicates the location of the 
threshold number no at JT = 0.02 (0.03). 



magnetic field {ip = 1^/2). In other words, the spectrum 
is flattened into the electric field plane (Fig. [3]). 

Eq. (|3.45[) features a factor \/Q va its exponent. Conse- 
quently, it is applicable literally only in the special case 
Q oi the tunneling regime. In this case the spectrum 
has a sharp peak at = (Fig. [3]), so that the positron 
and electrons are overwhelmingly emitted in the same 
direction. We expect though that the spectrum should 
exhibit similar main features within the general tunnel- 
ing regime (arbitrary C/, <C 1): a peak in the = 
direction and flattening into the electric field plane. 

The expression for the angular distribution in the cir- 
cular polarization case (|4.46p gives the mean polar angle 
of positron emission e = V2/^. Due to the ^V-^ factor 
in the exponent, this distribution is very sharp. That is 



the width of the distribution is 



56* - ^ « 1 . (6.8) 



Hence this distribution, which is symmetric in azimuthal 
angle t/j, is a thin- walled cone with axis of symmetry being 
the light wavevector and with conical angle -y/S/^ (Fig. 
S]). Note also that ip- — (p+ — tt at the momentum saddle 
point (see Eq. I4.15p . so the electron and positron are 
emitted in antipodal directions along this cone. 
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FIG. 3: Cross-sectional polar plots of the angular distribution 
(normalized to unity at its maximum) for the linear polariza- 
tion case for various fixed azimuthal angles ip. The polar 
coordinates are {dW/dQ,9), where the polar angle 9 is mea- 
sured from the direction of the light wavevector k, indicated 
by the arrow. Plots are shown for: k - £ or tp = plane (solid 
line); k - Ti or ip = n /2 plane (dashed line); and (p = n/A 
plane (dotted line). Here uj = 0.002m and = 0.01 so that 
^ = 5 and Q = 0.25. The distribution has a sharp peak at 
9 — 0, and is suppressed in the magnetic field direction. 



FIG. 4: Cross-sectional polar plots of the angular distribu- 
tion (normalized to unity at its maximum) for the circular 
polarization at arbitrary fixed azimuthal angle ip. The polar 
coordinates are {dW/dQ,9), where the polar angle 9 is mea- 
sured from the direction of the light wavevector fe, indicated 
by the arrow. Here lo — 0.002m and = 0.02 (solid line) 
or .7^ = 0.03 (dashed line) so that ^ = 10 or 15 respectively. 
The plots clearly exhibit very sharp peaks. Due to cylindrical 
symmetry the full distribution is obtained by rotating about 
k to produce a thin-walled cone. 



VII. CONCLUSION 

We considered the electron-positron pair production by 
Coulomb and laser fields in the tunneling regime describ- 
ing the process in simple analytical terms, and deriving 
expressions for angular and energy distributions. The en- 
ergy of the pair proves to be high, well above the thresh- 
old for the pair creation for any polarization. The angu- 
lar distributions strongly depend on the polarization. For 
the linear case both leptons move predominantly along 
the direction of the laser beam. For the circular polar- 
ization the pair is distributed in a thin-walled cone with 
axis of symmetry along the direction of the light beam, 
and with small conical angle. The leptons follow antipo- 
dal directions on the cone. The vicinal approximation 
suggested in this work proves to be very convenient for 
analytical calculations. 
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APPENDIX A: SPIN FACTORS 

In all polarization cases, the differential pair creation 
rate dW includes a spin factor 

K = J2\IC{to)\\ (Al) 

where IC{t) is defined in Eq. (|2.30p . and the subscript '0' 
indicates evaluation at the appropriate saddle points in 
the time and momenta variables. 

In the linear polarization case, from Eq. (|2.25p one 
finds that at the saddle points ti and Pq in the ^ ^ 1 
regime 

Q; = Q- = 1 + *^(7"-7')7' = Qo. (A2) 

It is straightforward to show that the matrix Qq has the 
property Qo ~ Qa- The spin factor is then 

^=^El^P-QoT°Qo"-P+l'' (A3) 

^ s± 

where at the saddle points ti and Pq, the lepton en- 
ergy So = 3m/2/-\/2 by Eq. (|2.16p . However, the free 
Dirac spinors u±p^ satisfy {jPo±)u±p^ = mu±p^ (see 
Eq. I2.13p ). Hence we may write the spin factor as a 
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trace of projection operators 

ilPo-) - mjQoj^Qo 



9 4m2 

X ((7P0+) +to)Qo7°<9o 



(A4) 



In the circular polarization case, once more from Eq. 
(|2.25p one finds that at the saddle points to and Pq in 
the ^ ^ 1 regime 



Ql = Q_ = 1 



(A5) 



This time, Qq ^ Qq, and the lepton energy Eq = 
I ^2. It follows then from Eq. (jA3p by similar rea- 
soning to the linear polarisation case that 



K 



2 1 



rTr 



((7P0-) - "i)Qo7°Qo 



^4 4m2 
X ((7P0+) + "i)Qo7"Qo 



2 



(A6) 
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